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Abstract 

The dynamics of perfect fluid spacetime geometries which exhibit Local Rotational Symmetry 
(LRS) are reformulated in the language of a 1 + 3 "threading" decomposition of the spacetime mani- 
fold, where covariant fluid and curvature variables are used. This approach presents a neat alternative 
to the orthonormal frame formalism. The dynamical equations reduce to a set of differential relations 
between purely scalar quantities. The consistency conditions are worked out in a transparent way. 
We discuss their various subcases in detail and focus in particular on models with higher symmetries 
within the class of expanding spatially inhomogeneous LRS models, via a consideration of functional 
dependencies between the dynamical variables. 



PACS number(s): 98.80.Hw, 04.20.-q, 04.20.Jb, 97.60.-s 



* e-mail: H.van.Elst@maths.qmw.ac.uk 
^e-mail: ELLIS@maths.uct.ac.za 



1 



1 INTRODUCTION 



2 



1 Introduction 

In studying the geometry and dynamics of relativistic cosmological models two main streams have oc- 
curred: one [1, 2] focusing on the spacetime (Killing) symmetries of these models, and one [3, 4] focusing 
on their covariant properties arising from a 1 + 3 "threading" decomposition of the spacetime manifold 
with respect to an invariantly defined normalised timelike congruence. It is commonplace to determine 
the covariant properties of models when the prime focus of analysis is imposed spacetime symmetries, 
but the converse case has not been systematically addressed. This paper is one of a series that aims to fill 
this gap, that is, to determine how spacetime symmetries appear when models have been characterised in 
terms of their covariant description. There are two main cases to consider here. The first is when there is 
no continuous isotropy group: there exist at most discrete isotropies in the cosmological model. That will 
be the subject of a later paper. The second is when they are Locally Rotattonally Symmetric (LRS) [5, 6], 
that is there exists a continuous isotropy group at each point, and consequently, a multiply-transitive 
isometry group acting on the spacetime manifold . Such models are the concern of this paper, which aims 
to covariantly characterise all cosmological models that are LRS and have a perfect fiuid matter source. 

LRS spacetime geometries have been discussed many times in the literature, usually in terms of a 
local coordinate or an orthonormal frame formulation. It is known that in the cosmological context, the 
isotropies around a point can occur as a 1-D or 3-D subgroup of the full group of isometries, which leaves 
the normalised 4- velocity u' / c of the matter fiuid fiow invariant. The latter case are the everywhere 
isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime geometries, that are well under- 
stood (and are the standard models of cosmology). The former are anisotropic and in general spatially 
inhomogeneous (the spatially homogeneous ones are characterised as subcases of the Bianchi models in 
[7, 8]). We will consider their nature and classification in detail in what follows. 

In the covariant 1 -|- 3 "threading" picture the invariantly defined normalised 4- velocity u' / c (such 
that Ui/cu' /c = —1) of the matter fiuid fiow and a tensor projecting orthogonal to it, defined by h'j := 
6'j + u' /cuj/c, are the key ingredients. We assume that the matter content of the cosmological model at 
hand, which in the fiuid description in general has stress-energy-momentum tensor 

Tij = p hij + TTij + 2 g(i/c Uj-j/c + /j, Ui/c Uj / c (1) 

when decomposed into its irreducible parts with respect to /c, has the characteristic features of a 
perfect fiuid (vanishing energy current density, = 0, and vanishing anisotropic pressure, tTjj = 0), 
and that the matter fiuid is isentropic, so that its total energy density jj, and its isotropic pressure p 
are related via a barotropic equation of state p = p{n) [9, 10] (with ^ ^ unless stated otherwise). 
Through this assumption a lot of relativistic thermodynamical processes of physical interest involving 
viscous effects and entropy production are excluded. Nevertheless, these simplifications on the side of 
the matter source provide the possibility of a neat treatment of the geometrical properties of the above 
mentioned cosmological models with some spatial isotropy. We also exclude the vacuum case, that is, we 
assume {ji-l-p) > 0. The usual definitions of the kinematic quantities (see [3, 4]) are employed and we 
follow the sign and index conventions of Kramer et al [1]. 

When the spacetime is LRS, there exists a unique preferred spatial direction at each point, covariantly 
defined for example by either a vorticity vector field, an eigendirection of a rate of shear tensor field, or 
a non-vanishing non-gravitational acceleration of the matter fiuid elements. This direction constitutes a 
local axis of symmetry: all observations are identical under rotations about it, and in particular they are 
the same in all (spatial) directions perpendicular to that direction [5, 6]. Thus there exists a preferred 
spacelike unit vector field e' 

a u'/c = , Cic} = 1 => e' (Vj-ei) = . (2) 

Due to the local rotational symmetry, all covariantly defined spacelike vector fields in the spacetime must 
be proportional to e' , that is, on defining as usual (?yi234 = —\/—9 ) 

ViUj/c := a(^ij-j/c+ ^(Q/c)hij - rnjkiui'' /cu' /c - Ui/ciij /c"^ , (3) 
with a\ = and ui* /c := — 1/2 rf^^^ ujj]./ cui/ c O i^ij /c = — rjij^iuj^ / cu^ / c, we must have 



u'/c'^ = {u/c'^)e\ h^iVjH = n' Ci , h^iVjp = p' Ci , 
cu'/c = {oj/c)e\ hhVj{Q/c) = {Q/c)'ei. 
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This applies equally to the derivatives of e': using (2), 

h){e^y := h)e' {Vke') =0, (5) 

h) {e^y/c := h) u^/c (Vfce^') = , (6) 

where a prime denotes the covariant derivative along the preferred direction , and a dot the covariant 
time derivative along the fluid flow lines. Hence, the preferred spacelike direction is geodesic in the 
local 3-spaces orthogonal to / c and Fermi-propagated along the matter fluid flow lines. The spatial 
rotation ("twist") of must be proportional to itself, that is 

ry'^'"(Vje,)«,/c= -fee', (7) 

where k denotes the magnitude of the spatial rotation of : k := \ rf^^^ jCk) ui/ c\. We deflne a as the 
magnitude of the spatial divergence of : 

a := hh (Vj-e') => Vie' = a + {ii/c^) . (8) 

A unique spacelike iracefree symmetric tensor field Cij is deflned from e' by 

which has the properties that eij =0, eij = €{, e\ = 0, e^j = 1/2 {eij + hij) and e^i = 3/2. 
We deflne a tensor projecting orthogonal to both and / c by := /ij — Cj. The spatial divergence 
and the covariant time derivative of Cij are 

h\ h^i (Vj-e") = f h^k {Vje') = f a e' , (10) 

h'ih'j{e,iy/c = 0. (11) 

Again, due to the assumed LRS symmetry of the spacetime manifold, a//tracefree symmetric tensor flelds 
have to be proportional to Cij. Thus, with the deflnitions 

Eij := Cikjiu''/cu'/c = E(ij-f , (12) 

Hij := -^Viki^C'"'jnu'/cu"/c = H(ij) , (13) 

which imply E'i = 0, H'i = 0, we have for the rate of shear tensor and the "electric" and "magnetic 
parts" of the tracefree Weyl conformal curvature tensor respectively 

where we deflne the squared magnitudes cr^ := 1/2 a) > 0, E'-' := 1/2 E) Eh > 0, i?^ — 1/2 H) H{ > 
0. Note that a, E and H as deflned in Eq. (14) can be either positive or negative. It follows that all 
solutions are of Petrov type D if the Weyl curvature tensor is non-zero [11, 1]. Thus we now flnd that 

ViUj/c = {(t/c) Cij + i (6/c) hij + (w/c) S[ij] - (m/c^) Ui/ccj , (15) 

with the antisymmetric tensor deflnition s^ij] := — rjijki / c. Further, on taking the orthogonality 
relation (Vjej) = into account, we can express the total spatial projection of the covariant derivative 
of as 

h\ h'j (Vfce,) = f pf^ijj + I S[ij] . (16) 
Expanding out the projection tensors and remembering (2), this leads to 

ViBj = f + I S[ij] + [ (a/c) + i (e/c) ] Ci Uj/c - {u/c"^) Ui/cuj/c . (17) 

From the point of view of the Cartan-Karlhede equivalence problem formalism for invariantly clas- 
sifying different spacetime geometries [12, 13], in the covariant 1 -|- 3 approach to LRS spacetimes we 
condense rotational isometries in the introduction of the spacelike unit vector held . To cover transla- 
tional isometries we can specify a maximal set of four generalised essential coordinates, which we then 
can investigate for the existence of functional dependencies [12]. For (non- vacuum) models with non- 
vanishing rate of expansion (0/c) we choose this set to be constituted by the four covariantly deflned 
scalars ^4 := { (6/c) , {<J / c) , E }. 

LRS perfect fluid spacetime geometries found frequent application in the literature in for example the 
modeling of the dynamical processes surrounding the formation of relativistic stars or galaxies ([14], [15], 
[9], [10]). Irrespective of the large degree of idealisation of the physics underlying their evolution, they 
have proved to be a valuable test ground for more complex astrophysical scenarios. 
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2 The Equations 

Using the relations introduced in the previous section, we can express all the covariant equations of Ref. 
[4] as equations for covariantly defined scalars /. These are in our case n, p, (0/c), (cr/c), (ui/c), (u/c^), 
a, k, E, and H . 

2.1 The Ricci identities 
2.1.1 Time derivative equations 

{Q/cy/c = -l(e/cf + {u/c'y + a{u/c') + {u/c'y-2ia/cf + 2iij/cf 



- ^ + 3p) (18) 

(co/cy/c = -l(e/c)(co/c) + ^(a/c)(co/c) + !l{u/c') (19) 

(a/cy/c = ^Au/c-)'-^^{u/c-) + ^^{u/c-f-^^{./cf-^^{a/cf 

-l{Q/c){a/c)-E (20) 

2.1.2 Constraint equations 

= (e/c)'-V3(<7/c)'-2^a(,7/c)-|fc(c^/c) (21) 

= (w/c)' + a(w/c) - (m/c^) (w/c) (22) 

H = -V3(M/c2)(w/c) + ^a(w/c) + |A;((T/c) (23) 

2.2 The Bianchi identities 

2.2.1 Time derivative equations 

E/c = -^{|l + p){(T/c) + ^/lE{(T/c)-{Q/c)E+lkH (24) 

H/c = V3H(a/c)-(e/c)H -^kE (25) 

11/ c = -{,j + p){e/c) (26) 

2.2.2 Constraint equations 

E' + laE = ^^' + 3H{cu/c) (27) 

H' + laH = -V3^{^ + p){iu/c)-3E{iu/c) (28) 

p' = -{^, + p){u/c') = 1^/ (29) 

2.3 Vanishing vorticity 



When the vorticity vanishes, (uj / c) = 0, for LRS perfect fiuid spacetime geometries the Gaufi equation 
[15], which describes the intrinsic curvature of spacelike 3-surfaces embedded orthogonal to a matter fiuid 
fiow into the spacetime manifold, can be reduced to 

'R^j = 7S [ ^ - i (-^/c) + ^ i^/c)' ] e., + i [ 4 ^ - i (e/cf + 2 (a/c)' ] h,, . (30) 

The trace of this equation is 

3R = 4^fi-l(e/c)' + 2(a/c)' , (31) 

and constitutes a generalised Friedmann equation. From this expression and the relations given by the 
Ricci and Bianchi identities one can also derive the conformally invariant 3-Cotton-York tensor defined 
by [16, 17] 

3C('J') := - /i^/^ r?'"™ Um/c ^Vk [ ^R^i - \ h^i ] ■ (32) 
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2.4 The equations for e' 

Given the decomposition (15), one obtains the Ricci identity for u' / c and so the equations listed above. 
Analogously, given (17), one can obtain the Ricci identity for e' and a corresponding set of equations. 

2.4.1 Ricci identities 

Now from the Ricci identities we can construct 

Vi {Vje') - {Vie')' - Rij e' = , (33) 

which gives 

= a' - I (e/cf - ^ (e/c) (a/c) + | (a /of + E + ^ - 1^ + I ^ , (34) 

and 

u^/cVi (Vj-e') - (Vie') '/c - Rij u'/ce^ = , (35) 

which gives 

= a/c + f (e/c) - ^ (a/c) - | (S/c) {u/c') + ^ (a/c) (u/c') - k (uj / c) . (36) 

2.4.2 Jacobi identities 

Next we consider the Jacobi identities or first Bianchi identities. Here we obtain from 

R[ijk]ie^ = ^ O V[iVjefc] = 0, (37) 

after contracting with jf^^^ ui/c, 

Q = k' + ak-l (e/c) (uj/c) - ^ (a/c) (uj / c) , (38) 

while contracting with rf^^^ ei, we get 

= A;/c+ iA;(e/c) - ■^A;((t/c) . (39) 

3 The Consistency Equations and their Different Cases 

We now turn to derive and investigate the consistency equations for the above set of covariant evolution 
and constraint equations for LRS perfect fiuid spacetime geometries. 

First we note that for any scalar invariant /, by the LRS symmetry, 

Vi/ = /'ei-//c«i/c, (40) 

because the spatial derivatives in directions perpendicular to must be zero. Taking the covariant 
derivative of this relation gives 

ViVJ = (Vjf) a + /' (VjCi) - {Vjf/c) Ui/c - f/c(VjUi/c) . (41) 

Multiplying by ry-*'*^' ej, m;/c, and using (uj/c) e' = — 1/2 ry'-**^' (Vj-m^/c) m;/c and (7), we find the important 
result 

2f/c(Lj/c) = f'k . (42) 
Applying this first to (lo / c) and using Eqs. (19) and (22), and then to k and using (39) and (38), we find 

(uj/c)D = Q = kD , (43) 

where 

D := I (e/c) (C./C) - ^ (a/c) (c./c) - ak . (44) 
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Now (43) implies D = 0, unless k = = (uj/c), but then D is zero in any case, from its definition. Hence 
we always have 

D = O ak =l(e/c)(co/c)-^(a/c)(co/c) . (45) 

Also applying (42) to p we find 

2|^(e/c)(c./c) = fc(«/c^) . (46) 

Now to check consistency, we take the covariant time derivative along the fiuid fiow lines of each of the 
constraint equations (21), (22), (23), (27), (28) and (29) and demand that they vanish. Within this 
procedure we make frequent use of the relation 

ify/c = {f/c)' + {u/c') f/c - ^ (a/c) f - i (e/c) /' , (47) 

where / denotes any arbitrary, covariantly defined scalar quantity, and re-substitute, where possible, for 
all spatial and temporal derivatives of the fiuid and curvature variables from the evolution and constraint 
equations (18) - (29). From Eq. (27), we obtain 

= Vi^{}i+p)k [uj/c) , (48) 

and since we are assuming {jJ. -\- p) > 0, we obtain from this 

A;(w/c) = 0. (49) 

Now we can multiply (42) first by k and then by (w/c) to attain 

f'k' = = f/c{cu/cf , (50) 

which is always valid. Similarly, from (46) we obtain 

k^u/c^)=0 = f^^(e/c)(^/cf . (51) 

Applying (50) first to k and then to (uj / c) we find that k' k = and (w/c)'/c(w/c) = 0, which implies 
that always 

k' = = (Lo/cy/c . (52) 
Putting the first of (51) and the second of (52) into Eq. (19) and then using (45) we find 

ak = 0. (53) 

The preferred spacelike unit vector field e' in LRS spacetime geometries is thus either (i) orthogonal to 
spacelike 2-surfaces (k = 0), and can be derived from a scalar potential, or (ii) divergence-free (a = 0), 
and can be derived from a vector potential, or (iii) covariantly constant (a = = k), when regarded in 
the local 3-spaces orthogonal to the matter fiuid fiow. Finally, putting the first of (52) in (38) and using 
(53), and employing the first of (51) and the second of (52) in the evolution equation (19), this shows 
that 

(e/c) (co/c) = = (a/c) (co/c) . (54) 

Using these results in the rest of the integrability conditions referred to above, the only non-trivial one 
is the propagation equation for the magnitude of the acceleration of the fiuid, 

= (u/c') -/c - [1^ (e/c)] ' - (1^ - i) (e/c) {u/c') + ^ (a/c) {u/c') , (55) 
which implies no new constraints (in checking these conditions, note that (54) implies ((jJ / c) -i^ (©/c) = 



d jj, 

0). 

Hence, the set of covariant evolution and constraint equations describing LRS perfect fiuid spacetime 
geometries with (jj, -\- p) > Q is consistent provided that the conditions 

k(uj/c) = Q, ak = Q = (Q/c)(uj/c) = (a/c)(uj/c) , f / c(uj / c) = Q = f k , (56) 

are true simultaneously. (The case f/c = (u/c^) is not included in this notation, as (u/c^) is not the 
time derivative of a covariantly defined scalar.) Then also the covariant time derivatives of the constraint 
equations (34) and (38) vanish identically. Assuming (ji-l-p) > 0, there are thus three cases that can 
occur [5, 6]: 
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• LRS class I: (w/c) 7^ ^ k = (e/c) = (a/c) = , //c = . 
e' is hypersurface orthogonal, u' / c is twisting. 

• LRS class II: A; = = (ui/c) . 

e' and u' / c both are hypersurface orthogonal. 

• LRS class III: k ^0 => (ui/c) = a= (ii/c^) = , /' = . 
e' is twisting, u' / c is hypersurface orthogonal. 

We discuss these in turn, in each case considering first the generic situation and then the subcases that 
occur. 

It remains to check the consistency of the algebraic expression (23) for H with the constraint equation 
(28). When this is done and the conditions (56) have been imposed, one obtains 

0= [V3(«/c2)' + V3(«/c2)'-^(/. + 3p)-2i?](c./c) . (57) 

This is only of interest when (w/c) 7^ (see the following section). 

4 (uj/c) ^ 0: Rotating Solutions (LRS Class I) 

When (w/c) 7^ 0, we immediately find that A; = = (6/c) = (c/c) , so the models within this LRS 
class of solutions can neither expand nor distort. Further, the consistency conditions (56) show that 
//c = for all covariantly defined scalars /, that is dct is a timelike Killing vector field and all spacetimes 
within this LRS class will also be stationary. Thus, in this case there exists a G4 multiply-transitive on 
timelike 3-surfaces. Non-zero quantities in general are n, p, (u/c^), a, E and H. 
The set of equations one needs to solve in this LRS class is given by 

(«/c')' = - a {u/ c^) - {u/ c^f -2{uj/ cf + ^{,i + ?,p) (58) 
{uj/c)' = -aiij/c)+{u/c')iij/c) (59) 
{u/c') = -|^/,'/(/,+p) (60) 

a' = -^ + a{u/c^)+2ii,/cf-2^i^ + p), (61) 

which follow from (18), (22), (29) and (34) respectively. The magnitudes of the "electric" and "magnetic 
parts" of the Weyl conformal curvature tensor in these models are given algebraically by a combination 
of Eqs. (20) and (18) and Eq. (23) respectively, which are 

E = -^a{u/c')-V3ic./cf + ^i^ + 3p) (62) 

H = -V3{u/c'^)(uj/c) + ^a(uj/c) . (63) 

With Eqs. (58) and (62) Eq. (57) is identically satisfied, as are the constraint equations (27) and (28) 
on substitution of (62) and (63). 

When the acceleration is non-zero, by combining Eq. (58) with (60) we can derive a coupled second- 
order ordinary differential equation for the spatial distribution of the energy density, which is 

1^ + + a 1^ - (2 1^ + 1) 1^ i, + p) 

-2{,i + p){uj/cf + ^{,i + p){,i + ?,p) = 0. (64) 



Algorithm: The boundary data, which can be specified freely on a timelike 3-surface orthogonal to e', 
are the values of ji, (w/c) and a. We assume a given equation of state p{n). {ii/c"^) is then given by Eq. 
(60), while E and H follow from Eqs. (62) and (63). 

For a treatment of this LRS class in terms of local coordinates and a metric tensor field refer to 
Stewart and Ellis [6], Eq. (2.8), and Kramer et al [1], Eq. (11.4). 
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4.1 Solutions with a = 

A subclass of solutions with a = exists. However, for consistency they demand an equation of state 
of the form = — 1/3// + const, which is usually dismissed as unphysical. One ordinary differential 
equation describing the spatial distribution of the total energy density /j, remains to be solved. It follows 
from Eq. (58) and is given by 

fi" - 1 V (/. + p) + 3 ^ (/.^ - p2) = . (65) 
One then obtains the following expressions for the remaining non-zero quantities: 



(66) 



Models of this kind can be regarded as (non-physical) generalisations of the Godel LRS case (see subsection 
4.5 below). 

4.2 Solutions with p = 

In the subclass of dust models, p = O (u/c^) = 0, deriving from Eq. (61) there is again one ordinary 
differential equation for jj, which remains to be solved: 

//"-|//'V/^-2^/^' = . (67) 

Then it follows that 



a = -i..' 

2nG 



E 



2tiG .. TT _ Vl /'27rG\l/2 , , 1/2 (^^) 



Models with these properties can also be regarded as generalisations of the Godel LRS case (see subsection 
4.5 below). 

If we were to impose the dynamical restriction (m/c^) = instead of p = 0, Eq. (60) would be solved 
by either ji' = 0, which leads to /' = (subsection 4.5 below), or = p = const, which leads to 
slight generalisation of Eq. (67). 



a 



4.3 Solutions with i7 = 

The dynamical restriction H = implies from Eq. (63) that a = 2 (u/c^). Then from Eqs. (58) and (61) 
we find that 



[uj/cf = -{ulc-'f + l^{ji + 2p) . (69) 



Finally, consistency of Eq. (59), on using Eq. (60), demands that = 1 p{^) = ^ -\- const. 
4.4 Solutions with = 

The dynamical restriction E = Q implies from Eq. (62) that 

(c./c)^ = -f («/c^)+|f (/. + 3p) . (70) 
Consistency with Eq. (59) then requires that 



2 ' 3a ci -t' ' 3a L 3 ^dp. 

~^"3c 

in order to obtain a purely "magnetic" solution. 



{u/c 



+ 3^(/" + 3p) = 0, (71) 
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4.5 /' = 0: Godel's rotating model of the Universe 

If we impose the additional condition that /' = 0, the symmetry group is a G5 multiply-transitive on 
the full spacetime manifold, which consequently is homogeneous. As the spatial derivatives of all non- 
zero scalar quantities / vanish, we have immediately from Eq. (60) that the matter moves geodesically: 
(m/c^) = 0. However, since no spacelike 3-surfaces of constant / exist, which are the group orbits of a 
simply-transitive G3, models of this LRS subclass are not spatially homogeneous (see Section 5 below). 
Equation (59) shows that a = 0, and from the algebraic equation (63) we obtain that the magnitude of 
the "magnetic part" of the Weyl conformal curvature tensor is zero, H = 0. All scalar quantities / are 
constant on the spacetime manifold, so the remaining equations are purely algebraic. From Eqs. (58), 
(61) and (62) we find 

2{u,/cf = ^{iJ + 3p) = 2^{i, + p) (72) 
E = -V3(c,/cf + ^(^ + 3p) . (73) 

It can easily be seen that this system of algebraic equations is consistent only, provided the equation of 
state is p(/u) = /U, that is for stiff matter. An equivalent configuration would be a combination of dust 
matter, p = 0, and a negative cosmological constant, A < [15], which gives Godel's rotating model of 
the Universe [18]. 

In summary, we observe that LRS class I of perfect fiuid spacetime geometries with equation of state 
p = pIji) contains a variety of stationary solutions, that, apart from the Godel model, are differentially 
rotating. Most of them, however, are of minor interest for astrophysical and cosmological purposes, as it 
is difficult to see how, given non-zero vorticity, the geometry of a star model could correspond to exact 
rotational symmetry about every point. 

5 k ^ 0: Homogeneous Orthogonal Models with Twist (LRS 
Class III) 

When k ^ 0, the consistency conditions (56) demand that /' = 0, (1^ / c) = = a: a// spatial derivatives 
vanish and it follows immediately from Eq. (29) that the matter in these models moves on timelike 
geodesies, that is (m/c^) = 0. Thus, since u^/c is normal and geodesic, all scalars / are spatially 
homogeneous and there exists a G4 of isometries multiply-transitive on spacelike 3-surfaces orthogonal 
to u'/c. That is, the spacetimes themselves are (orthogonally) spatially homogeneous (OSH) [7]. The 
non-zero quantities in the generic case are p, (0/c), (cr/c), E and H. Note that, apart from the 
class-defining spatial rotation of e' , all vectorial quantities are zero. From the 3-Ricci identities and the 
Gaufi equation for e' [15] the 3-Ricci curvature tensor of the spacelike 3-surfaces orthogonal to u' / c can 
be determined to be 

+ i [ I - 2 V3 + ^ (e/c) (a/c) - 2 (a/cf ] h,, , (74) 
the trace of which, when combined with Eq. (31), yields the generalised Friedmann equation 

^i? = ^^k"^ -2V?>E+-^{Q/c){a/c)-2{a/cf 

= i^^-lie/cf + 2ia/cf . (75) 

The magnitudes of the "electric" and "magnetic parts" of the Weyl conformal curvature tensor are 
determined algebraically from Eqs. (75) or (34) and Eq. (23) respectively as 

^ = i7^(0A)' + H0A)Wc)-;|Wc)' + f ^'-2^/" (76) 

H = Ik (a/c) , (77) 

and a time evolution equation for the dynamical variable k follows from Eq. (39). We emphasise the fact 
that throughout this LRS class H ^ 0, unless we deal with a spatially isotropic situation (see subsection 
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5.3 below). Using Eqs. (25) and (77) one can show that the conformal 3-Cotton-York tensor in this LRS 
class is given by (see also [19]) 



= - ^ h'/^ [H/c+lH (e/c) - ^ F (a/c) \ e,, . (78) 
The set of dynamical equations describing this class of expanding LRS models reads 



{Q/cy/c = -l(e/cf-2{a/cf-^{i, + 3p) (79) 

1,/c = -(j, + p)(e/c) (81) 
k/c = -lk{Q/c) + ^^k{a/c) . (82) 

We have tested for consistency of this set of equations with the evolution equation (24), which, with Eq. 
(76), is identically satisfied. An evolution equation for the 3-Ricci scalar, Eq. (75), can be derived: 

[3/?] ■/c=-|(e/c)3i?+^(^/c)3ij_V3fc2(^/c) . (83) 

Alternatively, in the face of our aim of using the set S^^ = {jj,, (6/c) , (c/c) , £" } as the (maximum set 
of) four generalised essential coordinates for the investigation of functional dependencies in models with 
non-zero rate of expansion, one can solve Eq. (34) for k'^ instead, obtaining 

fc2 = -f {Q/cf - ^ (e/c) {a/c) + I {a/cf + ^ + | . (84) 

As A;^ > is demanded, this leads to an algebraic restriction on the values of (0/c), (c/c) and E . The 
set of dynamical equations is then 

(e/c)yc = -i(e/c)^-2(^/c)^-^(/. + 3p) (85) 

Wc)yc = -^(^/c)^-|(e/c)(^/c)-i? (86) 
E/c = ^{^|l-p){(T/c)+A^/lE{(T/c)-{Q/c)E 

- {Q/cf {a/c) - V3 (e/c) {a/cf + 6 {a/cf (87) 

/i/c = - {^ + p){e/c) . (88) 

We have tested for consistency of this set of equations with the evolution equation (82), which, with Eq. 
(84), is identically satisfied. 

Algorithm: The initial data, which can be specified freely on a spacelike 3-surface orthogonal to the 
matter fiow, are the values of {Q/c), {a/c) and Given the equation of state p(//), all time derivatives 
are determined as well as k from Eq. (84) and subsequently H from Eq. (77). 

If local comoving coordinates are chosen, the line element can be cast into the following form (cf. Eq. 
(2.8) of Stewart and Ellis [6] and Eq. (11.4) of Kramer et al [1]): 

ds'^ = X'\ct) [ dx - h{y) dz f + Y'^{ct) [ dy'^ + Y^'^y) dz'^ ] - d{ctf . (89) 

Then one immediately obtains relations for the following scalars: 

, , X/c Y/c , , , I ( X/c Y/c\ , , 

(e/c) = -f + 2^ Wc) = ^(^-f -^j . (90) 

This enables one to integrate Eq. (82), which gives 

X 
Y 

where C\ denotes an integration constant. The solution to Eq. (83) is given by 



^ = Ci — , (91) 



3i 



>ij=24-^^, (92) 
y2 2 Y 
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where C2 denotes a further integration constant, and thus Eq. (75) constitutes a first integral. Conse- 
quently, one is left with only three essential dynamical equations, these being Eqs. (75), (80) and (81). 
Equivalent is the set of equations 

/ ■ \ 2 

-2^P (93) 
2 ^ 1^ (94) 




C2 


3C7i2 


X2 


Y2 


4 


y4 


+^ 






~ y2 


4 


y4 


Y/c2 




X2 


Y 




y4 



-2^p. (95) 

This set of equations was given in Ref. [5] for p = 0. Additionally we want to point out that, in contrast 
to C'l, no invariant geometric meaning can be associated with the constant €'2- In fact C2 is just a 
constant of integration, resulting from the fact that (94) is a first integral of (93) and (95); so any two of 
these equations implies the third. One obtains a form of the field equations independent of C2 from (95) 
and the difference between (93) and (94). 

Exact solutions to the dynamical equations within this LRS class have been discussed in Kramer 
et al, chapter 12.3 [1]. The different underlying simply-transitive subgroups G3 of the possible G4 
isometry groups can be of Bianchi Type-II, Type-VIII/III and Type-IX [7]. For Type-IX one can have 

> 0, which means can change its sign in this class of models. Type-II and Type-VIII/III have 
^R < throughout their entire evolution and are only distinguished by different values of the sectional 
3-curvature. Dynamically they are identical. 

5.1 Dust solutions 

Because the fiuid fiow is geodesic in any case, the dust subcases are simply distinguished by marginally 
simpler evolution equations (p = in Eqs. (81), (93) and (95)). 

5.2 Solutions with E = ("Pure magnetic") 

Imposing the dynamical restriction E = 0, Eq. (87) yields for (cr/c) ^ the algebraic condition 

(e/c)' = 4^ (5/, - p) - V3 (e/c) {ale) + 6 {a I cf , (96) 

which in order to satisfy Eq. (85) restricts the equation of state through an algebraic condition for 
The value of A;^, as determined by Eq. (84), is 

fc' = |^(/^+p), (97) 

which is clearly consistent with the condition A;^ > 0, and thus with (c/c) we have from Eq. (77) 
that H ^ 0. Solutions of this kind are of purely "magnetic" character as regards the Weyl conformal 
curvature tensor. An example is provided by the self-similar Bianchi Type-II OSH LRS solutions of 
Collins and Stewart [20], which have a linear barotropic equation of state of the form p(/u) = (7 — 1)/U, 
and one obtains E = for 7 = 6/5 (see also [21]). 



5.3 The FLRW subcase 

If we demand that for (cr/c) ^ the spacelike 3-surfaces orthogonal to m'/c be of constant curvature, 
that is from Eq. (74) that 

E=^[(e/c)-V3(a/c)](a/c) , (98) 

one obtains from Eq. (87) an algebraic expression for (0/c)^, which, when substituted in Eq. (84), 
gives = and thus violates the definition of this specific LRS class. Thus there do not exist any 
shearing solutions in this LRS class with isotropic ^Rij ; the ^R > - FLRW models are the only LRS 
models with k ^ and 3-spaces of constant curvature, as follows from Eq. (75). These are invariant 
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under a Gq of isometries multiply-transitive on spacelike 3-surfaces. In this case there exists a 3-D 
family of rotational symmetries rather than one. This is the exceptional case where the spacelike unit 
vector field is not uniquely defined, and there is no covariant feature that picks out a preferred spatial 
direction. Nevertheless, we can still find local coordinate or orthonormal frame bases as before. In the 
local comoving coordinates of Eq. (89) we have that Y = X . 

The Einstein static model is the special FLRW case when (e/c) = 0, and thus A;^ = 8/3 ^ > 
from Eq. (84). Then the spacetime manifold is invariant under a multiply-transitive G7 of isometries 
and therefore homogeneous. From Eq. (85) we find that for consistency an equation of state of the form 
p(yu) = — 1/3// = const is required, which can be interpreted as a dust model {p = 0) with a positive 
cosmological constant, A > 0. 

In summary, the general OSH LRS class III perfect fiuid spacetimes with equation of state p = piji) 
are characterised by the existence of a G4 isometry group multiply-transitive on spacelike 3-surfaces, with 
a simply-transitive subgroup G3 which belongs to one of the various Bianchi types. Special cases are 
the ^i? > - FLRW models. Of particular interest is the possibility of constructing simple cosmological 
models with purely "magnetic" Weyl conformal curvature tensor. This offers a step to studying the 
underlying physical mechanisms which could generate solutions of this peculiar kind. 

6 k = Q = {uj / c): The Inhomogeneous Orthogonal Family (LRS 
Class II) 

When k = Q = (uj/c), there exist 3-surfaces orthogonal to the fiuid fiow in which there acts a G3 multiply- 
transitive on spacelike 2-surfaces orthogonal to e' . These are the spherically symmetric solutions and 
their generalisations with plane and hyperbolic 2-spaces. All members within this class of expanding (or 
in the time-reversed case contracting), spatially inhomogeneous LRS models have vanishing "magnetic 
part" of the Weyl conformal curvature tensor, as follows directly from the constraint equation (23), that 
is 

H = . (99) 

The non-zero quantities in the generic case are /i, p, (0/c), (cr/c), (m/c^), a, and E. Depending on which 
of these dynamical variables might be zero, a broad variety of different special cases arises, which we will 
discuss in detail in this and the following section. 

From the 3-Ricci identities and the Gauf5 equation for e' [15] the 3-Ricci curvature tensor of the 
spacelike 3-surfaces orthogonal to u' / c can be determined to be 

%j = _ 1 [ a' 2 A' ] - i [ 2 a' I - 2 A' ] hij , (100) 

where K denotes the (constant) GauiJian curvature ^i? := 2 K of the 2-D spacelike group orbits orthogonal 
to e' and u' /c. A further constraint is given by the generalised Friedmann equation 

= - [ 2 a' I - 2 A ] 

= 4^i,-l(e/cf + 2(a/cf . (101) 
From this we obtain for a' the expression 

a' = -^a' + K-2^^ + ^ {Q/cf - {a / cf . (102) 
This relation can be combined with Eq. (34) to give a purely algebraic expression for E , which is 

Then the set of equations describing all LRS models within this class, the metric tensor of which can 
always be diagonalised [1], is 

{Q/cy/c = -l(e/cf+{u/c'y + a{u/c') + {u/c'y-2ia/cf-^i^ + 3p) (104) 
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- (e/c) ia/c) - ^ {a/cf - ^ - ^ a2 ^ 4 A- (105) 

ii/c = -(/, + p)(e/c) (106) 

a/c = -^(e/c) + ^(a/c) + l(e/c){u/c^)-^(a/c){u/c^) (107) 

K/c = -iK {e/c) + ^K{a/c) (108) 

Wc)' = j=^{e/cy-la{a/c) (109) 

{n/c') = -f^^^'/i^ + p) (110) 



a = --a 



, K_2iEG^^l (e/c)2 - (a/cf (111) 
A" = -a A (112) 

We have tested for consistency of this set of equations with the evolution equation (24), which, on using 
Eq. (103), is identically satisfied. The evolution equation (108) derives from demanding preservation in 
time along the matter fiow lines of the constraint (102), whereas expression (112) is a consequence of the 
constraint (27). The covariant time derivative of the constraint (112) vanishes identically. 

By application of Eq. (112) one can show that the conformal 3-Cotton-York tensor for all models 
within this LRS class vanishes (see also [19]), 

^Cij = , (113) 

that is, the spacelike 3-surfaces orthogonal to u' / c are conformally fiat. The evolution equation for the 
3-Ricci scalar (101) is 

[^R]-/c = -^[2{e/c) + V3{a/c)][^R + 2a{u/c^)] + 2V3[K-^]{a/c) 

- I [ (e/c) -V^ia/c) ] [ {u/c^y + {u/^f ] . (114) 

Alternatively, in the face of our aim of using the set S^^ = {jj,, (e/c) , (c/c) , A } as the (maximum set 
of) four generalised essential coordinates for the investigation of functional dependencies in models with 
non-zero rate of expansion, one can solve the combined Eqs. (34) and (102) for K instead, obtaining 

A- = |^/.+ ^- ^A-i[ (e/c) - V3 {a/c) f . (115) 

Then the set of relevant dynamical equations is 

(e/c)yc = -\{Q/cf +{u/^)' + a{u/^) + {u/^f -2{a/cf -^{p. + ?,p) (116) 

Wc)yc = ^^{u/c^y -^^{u/^) + ^^{u/c^f -^^{a/cf -l{Q/c){a/c)-E (117) 

A/c = -^{ii+p){(t/c) + VIe{(t/c)-{Q/c)E (118) 
ii/c = -(^+p)(e/c) (119) 
a/c = -|(e/c) + ^(^/c) + |(e/c)(«/c2)-^(^/c)(«/c^) (120) 

Wc)' = ^(e/c)'-|a(^/c) (121) 

E' = -laE+^i,' (122) 

(«/c^) = -f^^^'/i^+p) (123) 



We have tested for consistency of this set of equations with the evolution equation (108) and the con- 
straint equation (112), which, with Eq. (115), are identically satisfied. 

Algorithm: The initial data, which can be specified freely on a spacelike 3-surface orthogonal to the 
matter fiow, are the values of /j, and (Q/c), while {a/c), A and a can be given at a point and their 
spatial distribution is determined by Eqs. (121), (122) and (124) respectively. Then, given a choice of 
the equation of state p{n), (u/c^) follows from Eq. (123), all time derivatives are known and K is given 
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by Eq. (115). 



If local comoving coordinates are chosen, the line element can be cast into the following form (cf. Eq. 
(2.8) of Stewart and Ellis [6] and Eq. (13.2) of Kramer et al [1]): 



ds'^ = X'^ix, ct) dx'^ + Y'^{x, ct) [ dy'^ + Y.'\y) dz'^ ] - F-'\x, ct) d{ct f 
Then one immediately obtains relations for the following scalars: 



(e/c) 

{a/c) 



2 Y' 




X Y 










Y 


F (x/c 


Y/c 


V3 I X 


~ ~Y~ 



{u/c 



LEI. 

X F 



Eqs. (108) and (112) can then be integrated to give 

with Ci an integration constant. Next, from Eq. (101), we obtain the differential expression 



2 

"X2 



Y 



X^Y^ 

X Y 



y2 



(125) 

(126) 
(127) 

(128) 
(129) 

(130) 
(131) 



In integrating the field equations Eq. (101) is commonly used in place of Eq. (104). 

Exact solutions to the dynamical equations within this LRS class have been discussed by Kramer et 
al, chapter 14.2 [1], for the spherically symmetric case (A' > 0). 

By employing the identity (47), from the set of equations (116) - (124) we can derive the following 
useful relations, which show the non-linear growth of spatial inhomogeneities within this LRS class (cf. 
Refs. [22, 23]): 



[{e/c)' y/c 



- I [ 2 (e/c) + V3 (a/c) ],,'-(j, + p) (e/c)' 

- [ (e/c) + 2^3 (t/c) ] (e/c)' + {u/c')" + a {u/c')' + 3 {u/c') {u/c')' 
+ 6a{a/c)' -{u/c') [ i (e/c)^ 



(132) 



2 

3-/3 



2 47rG 47rG 

3 c" 



(e/c) (a/c) + f{a/c)' + ^^E 



{u/c') - {u/c')' 



(133) 



[(^/c)']yc 



- ^ [ (e/c) + 2V3{a/c) ] (e/c)' + ^ {u/c')" - ^ {u/c')' 
+ V3 {ii/c') {ii/c'Y + 2^3 a (cr/c)^ + | a (e/c) (a/c) + | a A 

- 75 [ i {e/c)' + ^ (e/c) {a/c) + i {a/c)' + ^ A 



^ - i ^ - ^ (/. + P) (1^) " + f {u/c') - {u/c') ' 



E'\/c 
[a'\/c 



'-[2{e/c) + V3{a/c)]^l,' 



'-{^^ + p){e/c)' 



+ l^{fi+p)a{a/c) + [2a-{u/c')][{e/c)-V3{a/c)]E 
I [ (e/c) - V^{a/c) ] {u/c')' - ^ (e/c)' - ^ (e/c)' {a/c) 
+ I (e/c) {a/c)' + ^ {a/c)^ + ^ [ (e/c) - VI {a/c) ] 



(134) 
(135) 
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+ ^[2(e/c) + V3(./c)] [i? + 2^/.] 

- {u/c^) [ f (e/c) - ^ a (a/c) - | (e/c) (u/c^) + ^ (^/c) (n/c^) ] (136) 
['R']7c = -4[(e/c) + V3(l + |^)(^/c)] 

+ [ f (Q/^)' + ^3 (Q/^) (-^Z^) - f (-^Z^)' - 2 « -^^E-I^i,] (e/c)' 

+ 2a((T/c) [ (e/c) (a/c) + bVs (a/c)^ + 6 E ^ 

- I [ (e/c) -V^(a/c) ] {u/c^y - I a [ (e/c) + 2^3 (^/c) ] {u/c')' 

- 4 [ (e/c) - V^(a/c) ] («/c^) («/c')' - in/c') [ - ^ (e/c)' + ^ (e/c)^ (^/c) 

- I (e/c) (^/c)^ + ^ (^/c)^ - ^ [ (e/c) - V3 (^/c) ] 

+ I [ (e/c) + 2 V3 (a/c) ] a («/c^) + | ^ (7/^ + 9p) (^/c) + f ^ (e/c) 

- I (e/c) + 2^3 (^/c) ] + I [ (e/c) - V^(a/c) ] («/c^)' ] • (137) 

In scenarios of the formation of structure on cosmological distance scales it is well-established practice to 
study the evolution of Uneansed perturbations around an assumed FLRW background spacetime geom- 
etry. In their FLRW-linearised form the set of equations (132) - (137) simplifies significantly and could 
be used for example to investigate covariant and gauge-invariant spherically symmetric perturbations of 
an FLRW spacetime [22, 23]. The description of scalar (total energy density) perturbations, which are 
the only ones contributing at linear order, would be entirely covered by only two equations from this set, 
for example Eqs. (132) and (133), and the equations underlying the evolution of the FLRW background 
spacetime geometry (see subsection 7.2.1 below). 

6.1 Spatially inhomogeneous LRS dust models 

With the condition p = O (u/ c^) = imposed, the cosmological models within this LRS subclass are the 
Lemaitre-Tolman-Bondi spherically symmetric solutions [24, 25, 26] and their generalisations to spacelike 
2-surfaces with vanishing or negative Gaufiian curvature scalar A' [5]. These spatially inhomogeneous dust 
models belong to the so-called "silent" class recently discussed by Bruni et al [27]. The relevant equations, 
which we obtain from specialisation of the set (116) - (124), are 



(e/c)yc = 


-l(e/cf-2(a/cf-^ 




(138) 


(^/c)yc = 


-^^(a/cf-l(e/c)(a/c) 


- E 


(139) 


E/c = 


-i^^(a/c) + VlE(a/c)- 


-(e/c)E 


(140) 


/i/c = 


-li(e/c) 




(141) 


a/c = 


- 1 + 73 (-^/c) 




(142) 


(a/c)' = 


^,(e/c)' -la(a/c) 




(143) 


E' = 


3 „ J? , A-kG ..I 
2 « ^ + V3c^ 




(144) 


a' = 




- 3 I'^/C) - 75^- — - 3— 


(145) 



Algorithm: The tmttal data, which can be specified freely on a spacelike 3-surface orthogonal to the 
matter fiow, are the values of jj, and (e/c), while (a/c), E and a can be given at a point and their spatial 
distribution is determined by Eqs. (143), (144) and (145) respectively. Then all time derivatives are 
known and K is given by Eq. (115). 

Exact solutions to the dynamical equations within this LRS subclass have been discussed by Lemaitre 
[24], Tolman [25], Bondi [26], and Kramer et al, chapter 13.5 [1], for the spherically symmetric case 
(A' > 0). Ellis [5] discussed the more general cases. For dust, in the local comoving coordinates of Eq. 
(125), we have from Eq. (129) that A = 1. The set of equations (132) - (137) describing the non-linear 
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growth of spatial inhomogeneities simplifies considerably and reads: 

[,,'\/c = -I [2(e/c) + V3(^/c)]/.'-/.(e/c)' (146) 
[(e/c)']yc = -42^/,'-[(e/c) + 2V3(,7/c)] (e/c)' + 6a(,7/c)' (147) 
[(-^Z^)']'/^ = -^/^'-;fe[(e/c) + 2V3(^/c)](e/c)' + 2V3a(^/cf 

+ I a (e/c) ((t/c) + I a (148) 

[E'\lc = -l[2ie/c) + V^ia/c)]^^'-^^ie/cy + l^^aia/c) 

+ 2aE [(e/c)-V3(a/c)] (149) 
[a']-/c = - ^ (e/cf - ^ (e/cf (a/c) + | (O/c) (a/cf + ^ (a/cf 

+ [(e/c)-V3(^/c)] + ^[2(e/c) + V3(^/c)] [i? + 2^/.] (150) 

['R']-/c = -4[(e/c) + V3(a/c)]^i,' 

+ 2a((T/c) [ (e/c)((T/c) + 5^3 ((T/c)V6i; 1 . (151) 



The linearised set of equations corresponding to [22] easily follows. 
6.2 (c/c) = 0: The shearfree subcase 

An interesting subcase within LRS class II arises, if we demand that (c/c) = 0. Then we immediately 
obtain from Eq. (117) that the magnitude of the "electric part" of the conformal Weyl curvature tensor 
is given by 

^=75(«/^')'-i7i(«/^') + 73(«/c^)' • (152) 
Hence, in order to obtain a solution with E ^ 0, the fluid acceleration has to be non-zero, and we also 
assume a non-zero rate of expansion of the matter fluid. (The case with (0/c) = will be treated in 
subsection 7.1 below.) Next, from the constraint (121) it follows that 

(e/c)' = , (153) 

that is, the spatial distribution of the rate of expansion is homogeneous and thus constant on the spacelike 
3-surfaces orthogonal to u'/c. With Eqs. (152) and (110) the conditions deriving from Eqs. (133) and 
(134) for (c/c) = are equivalent: they provide the expression 

+ iu/c') [ i (e/cf + ^-l^,-la iu/c') - i (u/c^y ] , (154) 

which then ensures that the constraint (122) is solved identically. Using relation (110), Eq. (154) is 
a constraint equation for the spatial distribution of the total energy density /j, of third order, which is 
coupled to the constraint equation (124), giving the spatial distribution of the spatial divergence a. 
In the local comoving coordinates of Eq. (125) we have that 

(a/c) = O^F^ = F^^^ (e/c) = 3F^. (155) 
Then we can immediately integrate Eq. (132) to find 

= || , (156) 

where C2 denotes an integration constant. Using this result and Eqs. (110) and (129), an expression for 
4^ is established. 

Exact solutions to the dynamical equations within this LRS subclass have been discussed in Kramer 
et al, chapter 14.2.3 [1], for the spherically symmetric case (A' > 0) and p = p(/u). 
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It would be interesting to investigate, whether a solution with (0/c)'/c = (Q/c) = const, 
(ji + 3p) > and negative deceleration parameter q exists within this subclass of LRS class II, which 
is consistent. A solution of this kind would provide a very special spatially inhomogeneous inflationary 
cosmological model. 

6.3 a = 0: The non-diverging subcase 

If a = 0, the normals to the spacelike 2-surfaces spanned by the isometry group are non-diverging. The 
conditions for this dynamical restriction to be consistent are the following: for a = Eq. (120) demands 
that 

0=[(e/c)-V3(a/c)] {u/c') , (157) 

while Eq. (124) gives E algebraically. If (i) (Q/c) = \/3(a/c), inserting E into the constraint (122) 
requires that = 0, while Eqs. (117) and (118) give (jJ- + p) = and (// + p) (0/c) = respectively. 
Consequently, this case can be discarded. If on the other hand (ii) (m/c^) = 0, for f ^ we need to 
have p = 0. Then from Eq. (144) we get the condition 

(Q/c)' = 3 [ (e/c) - V^{a/c) ]-i ^ ij' , (158) 

and the covariant time derivatives along the matter fluid flow lines of the constraints (143), (144) and 
(145) vanish identically. That is, only dust provides an appropriate matter source for spatially inhomo- 
geneous models of LRS class II with non-diverging isotropy generator. These models are thus a further 
specialisation of the "silent" class [27] discussed in subsection 6.1. 

In conclusion, the generically spatially inhomogeneous LRS class II constitutes the largest class of 
solutions of perfect fluid spacetime geometries with equation of state p = p(/u). The models are in general 
time-dependent, and in the past they have been prime candidates for the theoretical and numerical 
description of star and galaxy formation processes as well as supernova explosions (see Refs. [9], [10], [14] 
and [15]). Due to their highly idealised spacetime symmetry properties the relevant dynamical equations 
become relatively simple and easily tractable. However, as outlined in the introduction, an isentropic 
matter fluid flow, resulting from a barotropic equation of state, will in general be too restrictive to 
realistically model for example explosions of stars during the late stages of their evolution. 

In the following section we discuss those subcases of cosmological models within LRS class II, which 
contain an additional (timelike or spacelike) Killing symmetry, apart from the multiply-transitive G3 
isometry group. 

7 Higher Symmetry Subcases: Hypersurface Homogeneous 
Models in LRS class II 

Given generic initial data in the spherically symmetric family, we obtain general spherically symmetric 
inhomogeneous solutions, and their plane symmetric and hyperbolic generalisations, as discussed in the 
last section. Now we consider those cases where the initial conditions chosen result in higher symmetry 
solutions. By this we mean in particular that through such a special choice a further Killing spacetime 
symmetry shall arise. As a preliminary we note that if we have for example spatial homogeneity, then 
necessarily any covariantly deflned scalar / is functionally dependent on say the total energy density /j, 
[28, 29]. That is 

/ = /(/«) (159) 

and so 

V,/ = ^ V,/. , (160) 

which implies the two relations 



i^,:=/.^-,V,/ = ^/.^-,V,,:=^X,, 



(161) 
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and 

Similar relations will hold for each pair of scalars, like for example { (0/c) , /}. Before discussing the 
spatially homogeneous subcases within LRS class II in detail, we first want to turn our attention to the 
situation where the special choice of initial data gives rise to timelike 3-surfaces of homogeneity. 

7.1 //c = 0: The static subcase 

Imposing the geometrical condition //c = on the set of equations (116) - (124), that is, allowing for the 
additional timelike Killing vector field dct such that the group of isometries is a G4 multiply-transitive 
on timelike 3-surfaces, for (ji + p) > 0, E ^ 0, we immediately get from Eqs. (119) and (118) that 
(0/c) = = (cr/c). Then by combining Eqs. (116) and (117) we obtain for the magnitude E 

which, when substituted into Eq. (115), yields for the (constant) GauiJian curvature of the spacelike 2-D 
symmetry orbits of the G3 subgroup 

K = ^ + a{u/c^)-2^p. (164) 

In principle K can be positive, negative, or zero. Finally, the remaining equations of the set (116) - (124) 
lead to the set 

(u/c')' = -a{u/c')-{u/c'y + ^i^ + 3p) (165) 
(«/c2) = -|£/,7(/, + p) = -p'/(^^+p) (166) 

a' = -^ + a{u/c')-2^i^+p) . (167) 

With Eq. (163) the constraint equation (122) is identically satisfied. 

In the local comoving coordinates of Eq. (125) we have that X = X(x), Y = Y(x), F = F(x). In the 
case of modeling static, spherically symmetric matter configurations (where K > 0), however, it is more 
convenient to use Schwarzschild coordinates such that the line element assumes the form [15] 

ds^ = e^('') dr^ + {dd"^ sin^ d dip^) - e''^'''' d(ctf , (168) 

and the fiuid 4-velocity is given by / c = e'"^^ Upon combination with the Einstein field equations 
this leads to [15] 

(u/c^)=[l-'-^Y' ^[m(r) + ^r^p] , (169) 

where m(r) := JJ" /^{x) x^ dx. Inserting (169) into Eq. (166) then yields the Tolman-Oppenheimer- 
Volkoff structure equation, which describes the interior spacetime geometry of static, spherically symmet- 
ric relativistic stars [9, 10]. If we use expression (169) in Eqs. (165) and (167), we have to solve a coupled 
set of two ordinary differential equations for the two functions /u(r) and a(r), provided an equation of 
state p = pip) has been specified. 

In the remaining two subsections we want to address the spatially homogeneous subcases within LRS 
class II, which generally occur in two different ways. 

7.2 Orthogonal spatially homogeneous models (OSH) 

This is the higher symmetry family in which there is an extra symmetry acting on spacelike 3-surfaces 
orthogonal to u' /c, and with shear. Consequently, additional to (uj / c) = 0, we have /' = for all scalars, 
which necessarily implies from Eq. (123) that (u/c^) = 0, and also from Eqs. (121) and (122) that a = 0. 
Thus with /' = 0, Eq. (161) is identically satisfied because we have that 



Xi = Q , Fi = 



(170) 
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Equation (162) will follow from the dynamics of the system. It does not tell us anything about spatial 
homogeneity. 

The 3-Ricci tensor for OSH models of this subclass is given from Eq. (100) by 

^Rij = I K (hij - Cij) = Kpij , (171) 
while the generalised Friedmann equation is 

^R=2K = i^^-lie/cf + 2ia/cf . (172) 
The evolution of ^R is governed by the equation 

[3/?] ■/c=-i [2(e/c) + V3((T/c)]3i? + 2V3A'(,T/c) , (173) 
and from Eq. (124) it follows that 

The dynamical equations (116) - (124) reduce to the set 

(e/cy/c = -l(e/cf-2ia/cf-^i^ + 3p) (175) 

1,/c = - (i, + p)(e/c) . (177) 

The isometry group underlying OSH models within this subclass can be of Bianchi Type-I/VIIo (where 
the group orbits have isotropic 3-Ricci curvature with ^R = = 2 A'), or of a special case of Bianchi 
Type-III (with K < 0) [7], or instead it can be the Kantowski-Sachs (KS) case [30] (with K > 0), which 
is the unique case where the group of isometries G4 multiply-transitive on spacelike 3-surfaces orthogonal 
to m'/c does not contain a simply-transitive G3 subgroup [4, 31]. 

In the local comoving coordinates of Eq. (125) we have that X = X(ct), Y = Y(ct), F = I, and the 
solution to Eq. (173) is given by 

'^ = 'Y-.+^' (178) 

where C'l and C2 are integration constants. Thus Eq. (172) constitutes a first integral and is generally 
used in place of Eq. (175). 



7.2.1 FLRW 

Here we have that (m/c^) = O (c/c) = O A = 0. This is a subcase of the previous one, but it 
no longer follows that necessarily a = 0, a' = 0, as this is the exceptional case where the spacelike unit 
vector field e' is not uniquely defined. As mentioned in subsection 5.3 there is no covariant feature that 
picks out a preferred spatial direction. 

With the condition a' = — 2 K required to obtain spacelike 3-surfaces orthogonal to u' /c, which have 
isotropic curvature (cf. Eq. (100)), the Friedmann equation is given by 



3ij=6 [a'-^] =4^//-|(e/c)' . (179) 
The evolution and constraint equations for the non-zero variables are 

{e/cy/c = -i(e/cf-^i^ + 3p) (iso) 

/i/c = -(^ + p)(e/c) (181) 

A/c = -|A(e/c) (182) 

a/c = -f(e/c) (183) 

[3i?]yc = -|(e/c)3i? (184) 

A' = -a A . (185) 
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^R=^, (186) 



However, once the sign of was fixed, where either ^R > 0, ^R = or ^R < 0, only two of these 
equations are essential, which are commonly chosen to be Eq. (179) and Eq. (181). In the local 
comoving coordinates of Eq. (125) we have that X = X(ct), Y = X(ct) Z(x), F = I, and Eq. (184) can 
be integrated to give 

_C 

X 

where C is an integration constant. 

The Einstein static model is the special FLRW case when (0/c) = and thus 1/6% = K — = 
2/3^^// > 0. From Eq. (180) we then find that for consistency an equation of state of the form 
pIji) = — 1/3// = const is required, which can be interpreted as a dust model with a positive cosmological 
constant, A > 0, as was already discussed at the end of subsection 5.3. 

7.2.2 The dust subcases 

Because the fiuid fiow is geodesic in any case, the dust subcases are simply distinguished by marginally 
simpler evolution equations: p = in Eqs. (175) and (177). 

7.3 Tilted spatially homogeneous LRS dust models 

Here there exists an extra symmetry, but not in the 3-surfaces orthogonal to u' /c. Instead, the 3-surfaces 
of homogeneity are iiHed with respect to u' /c, that is u' / c does not provide the unit normal to these 
3-surfaces [8]. In the tilted case, with the cosmological conditions (0/c) 7^ 0, (// + p) > imposed (which 
we will assume), we will have Xi ^ 0, Fi ^ 0. Then equation (161) implies that their spatial directions 
are parallel: 

1^'^" XjFkUi/c={) . (187) 

However from the LRS condition we have Xi = ji' e,- and Fi = /' e,-, so (187) will be trivially satisfied. 
Thus we also need the magnitude information. We don't know fiji) (see (159), but we do know it must be 
the same in both equations, so we get a non-trivial magnitude relation by eliminating df(pi)/diJ, between 
the two equations (161) and (162), leading to 

l^/c Fi = f/c Xi O 12/c f = f/c n' . (188) 

which is another way of saying that the surfaces of constant jj, and / must be the same. The corresponding 
set of equations to Eqs. (188) will hold for every pair of scalar invariants /, defined from the kinematic 
quantities, matter variables, Weyl curvature tensor, and all their covariant derivatives. 

The key question in each case, related both to the Postlate of Uniformal Thermal Histories (PUTH) 
[32] , and to the Equivalence Problem series of papers [12, 13] , is how few such relations we have to impose 
before consistency conditions imply all the rest must be satisfied; and which variables should we choose for 
these relations. The weak conjecture is that 84^ = {jJ,, (6/c) , (c/c) , E } will suffice. The strong conjecture 
is that S2 = {jJ., I c) } will suffice. However the latter may be disproved by the PUTH counter-examples. 

It is known that the perfect fiuid cases of this class are of Bianchi Type-V (see King and Ellis [8]), 
studied in some detail in [33] and [34]. However for reasons of simplicity, we here focus on the situation 
where the equation of state of the matter is that of dust, p = O (m/c^) = 0. The exact solution 
to the field equations for this non-rotating tilted LRS SH dust model was given by Farnsworth [35]. 
It belongs to the "silent" class [27] of subsection 6.1. We assume the functional dependency relations 
(0/c) = (0/c)[//], (c/c) = (c/c)[//], E = E[ijl\, a = a[//] and K = A'[//], and that all of these dynamical 
variables are non-zero. Then, using the set of dynamical equations (138) - (145), from Eq. (188) the 
following five consistency equations need to be satisfied: 
First, 

ii/c (e/c)' - {Q/cy/c fi' = , (189) 



which can be solved to give 



. , ^// + i(e/c)' + 2(,T/c)' , , , 

(e/c)' = - -;(e/,) / . (190) 
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Second, 

which, on using Eq. (190), gives 



/«/c (cr/c)' - (a/cy/c = , (191) 



,' = -1 a,(e/c)(a/c) 

' ^-^/^-i7^[(e/c)-V3(./c)p 



Third, 



ii/c E' - E/c n' = , (193) 
which with Eq. (192) leads to the algebraic relation 

Fourth, the condition 

ii/c a' - a/c = , (195) 

which, after inserting from Eqs. (190), (192) and (194), leads to a higher-order multivariate polynomial 
equation in the variables (0/c), (cr/c), E and a, which is quadratic in a, and finally the relation 

/«/c A" - k/c n' = Q , (196) 

which, after inserting from Eqs. (190), (192) and (194), yields a second higher-order multivariate poly- 
nomial equation in the variables (0/c), (cr/c), E, that however does not contain the variable a. 

The next step to take is to consider the covariant time derivatives along the matter fiuid fiow lines 
of the constraints (189), (191), (193), (195) and (196), and work out the consequences which establish 
that these be preserved in time. Those of the relations for the pairs {ji, (Q/c) }, Eq. (189), (cr/c) }, 
Eq. (191), and {//,£"}, (193), vanish identically, when Eqs. (190), (192) and (194) are employed, while 
the relations for the pairs {//,a}, Eq. (195), and {//,A'}, Eq. (196), give two further higher-order 
multivariate polynomial equations, the first again quadratic in a, while the second again does not contain 
a. Combining them with the two obtained from the functional dependency relations (195) and (196), we 
have to solve a set of four simultaneous non-linear algebraic equations in the variables (0/c), (cr/c), E 
and a. We tackle this problem by means of application of the GROBIER package of the algebraic computing 
system REDUCE. The solutions we find in this way are 



(0/c) = , 


(a/c) -- 


= , 


(197) 




(a/c) = 




(198) 




(a/c) = 




(199) 




(a/c) = 


75 (Q/^) ' 


(200) 


^ = -576(6/^)' ' 


(a/c) = 


-76(9/^) • 


(201) 



However, all of them have to be discarded, as they lead to consequences for /j, and which violate 
the initial assumptions. The failure of solution (197) is obvious, while (198) with (194) gives a zero 
denominator in Eq. (192), and solutions (199) - (201) lead to = from Eq. (194). 

Thus, as one can show that the alternative assumptions a = 0, K ^ and a = 0, A' = also lead to 
contradictions, we conclude that for consistency Eqs. (189) - (196) have to be solved for A' = 0, a 7^ 
instead (see also [31]). That is, from Eq. (103) we have 

^=7ll'" + ^«'-i7^[ (Q/^) - ^ ]' ' (202) 

which simplifies our considerations significantly. With (202) the evolution and constraint equations (140) 
and (144) are identically satisfied. Now Eq. (190) is still valid, whereas from Eq. (191) we get 

3V^a^(e/c)(a/c) ^203) 



(0/c)-V3(,7/c)]2-|a2 
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which shows that for a = we necessarily have yu' = 0, and then successively /' = for all scalar 
invariants /. With Eqs. (190) and (203) inserted, Eq. (195) yields the algebraic expression 

— fi-e (e/c) -3{a/c) -4a — — 2 (204) 

L J [ (e/c) - V3(<t/c) ] -|a2 

while Eq. (193) is identically solved. In principle, we can invert this (quartic) algebraic equation to 
determine the dynamical variable a in terms of (0/c) and (cr/c). Note that for a = Eq. (204) 
reduces to Eq. (172) with K = 0. Note also that for = 4/3 [ (e/c) - V3(a/c) f, from Eq. (204) 
there occurs a singularity in the Ricci curvature tensor, and from Eq. (202) simultaneously in the Weyl 
curvature tensor. Equation (204) is compatible with Eq. (203), when Eqs. (143), (145), (190) and (202) 
are employed. 

It remains to investigate the covariant time derivatives along the fluid flow lines of the constraints 
(189), (191), (193) and (195). We find that they vanish identically, after the expressions (190), (203) and 
(204) have been substituted, and thus consistency of our assumptions and results for the tilted SH dust 
subcase within the generically spatially inhomogeneous LRS class II has been established. 

Algorithm: In the tilted SH dust subcase of LRS class II the tmttal data, which can be specified freely 
on a spacelike 3-surface orthogonal to the matter fiow, are the values of ji, (Q/c) and (cr/c) at a point. 
Then the values of a and E at that point follow from Eqs. (204) and (202) respectively. The spatial 
derivatives on the initial spacelike 3-surface orthogonal to u' / c of the non-zero dynamical variables are 
determined from Eqs. (203), (190), (143), (145) and (144), while their time derivatives come from Eqs. 
(138) - (142). K = 0, as demonstrated. 

Applying the steps outlined in this subsection to the case with non-vanishing pressure and hence 
non-zero acceleration should lead to similar results. We leave this investigation open for future work. 



8 Conclusion 

All three distinct classes of spacetime geometries that arise for perfect fiuids with equation of state 
p = pIji) and LRS spacetime symmetries as well as their different subcases, defined through further 
either dynamical or geometrical restriction, have been discussed in detail in the relevant paragraphs of 
the last four sections. A short description of their various fields of application was given at the end of 
each section. Here we present a final compact survey of the different cases we have been dealing with in 
this paper and include references to the (sub)sections in which their discussion can be found. 

This may be contrasted with the Table given at the end of [4], which characterises the subsets in a 
different way. 
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Table 1: Survey of the existing LRS perfect fluid spacetime geometries with barotropic equation of state 
p = pIji) [ (jj, + p) > ]. References to subsections which discuss special subcases are given. 





LRS class I 


LRS class II 


LRS class III 




^0 








k 








^0 


(e/c) 





•^0 


^0 






• 0, static (7.1) 




{a/c) 





•^0 








• (6.2) 


• 0, ^i? > - FLRW (5.3) 


(«/c2) 


•^0 


•^0 







• (4.2) 


• 0, dust/"silent" (6.1) 




a 


•^0 


• 7^ (includes tilted SH (7.3)) 







• 0, "unphysical" (4.1) 


• 0, "silent" or OSH/KS (6.3/7.2) 




E 


•^0 


•^0 






• 0, "magnetic" (4.4) 


• 0, FLRW (7.2.1) 


• 0, "magnetic" (5.2) 


H 


•^0 









• (4.3) 




• 0, ^i? > - FLRW (5.3) 



